Optimal thermal refrigerator. 
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We study a refrigerator model which consists of two n-level systems interacting via a pulsed 
external field. Each system couples to its own thermal bath at temperatures Th and Tc, respectively 
{9 = Tc/Th < 1). The refrigerator functions in two steps: thermally isolated interaction between 
the systems driven by the external field and isothermal relaxation back to equilibrium. There is 
a complementarity between the power of heat transfer from the cold bath and the efficiency: the 
latter nullifies when the former is maximized and vice versa. A reasonable compromise is achieved 
by optimizing over the inter-system interaction and intra-system energy levels the product of the 
heat-power and efficiency. The efficiency is then found to be bounded from below by C^ca = ^^_g ^ 1 
(an analogue of Curzon-Ahlborn efficiency for refrigerators), besides being bound from above by the 
Carnot efficiency ("c = — 1- The lower bound is reached in the equilibrium limit 6 —f 1, while 
the Carnot bound is reached (for a finite power and a finite amount of heat transferred per cycle) in 
the macroscopic limit Inn S> 1. The efficiency is exactly equal to Cca, when the above optimization 
is constrained by assuming homogeneous energy spectra for both systems. 

PACS numbers: 05.70.Ln, 05.30.-d, 07.20.Mc, 84.60.-h 



I 

o 



> 
m 

00 

in 
o 

ON 

o 



X 



Thermodynamics studies principal limitations imposed 
on the performance of thermal machines, be they macro- 
scopic heat engines or refrigerators [1], or small devices 
in nanophysics [2] and biology [3]. Let us recall three 
basic definitions applicable to any thermal machine tak- 
ing as an example a refigerator driven by a source of 
work: i) heat Qc transferred per cycle of operation from 
a cold body at temperature Tc to a hot body at tem- 
perature Th {Th > Tc). ii) Power, which is Qc divided 
over the cycle duration r. Hi) Efficiency (or performance 
coefficient) C, = Qc/W, which quantifies the useful out- 
put Qc over the work W spent by the work-source for 
making the cycle. The second law imposes the Carnot 
bound C ^ Cc = Tc/ {Th — Tc) on the efficiency of re- 
frigeration [1]. Within the usual thermodynamics the 
Carnot bound is reached only for a reversible, i.e., an in- 
finitely slow process, which means it is reached at zero 
power [1] . The practical value of the Carnot bound is fre- 
quently questioned on this ground. The drawback of zero 
power is partially cured within finite-time thermodynam- 
ics (FTT), which is still based on the quasi-equilibrium 
concepts [4] . For heat-engines FTT gives an upper bound 
V < VCA = I - \/Tc/Th (Curzon-Ahlborn, or CA ef- 
ficiency) for the efficiency rj at the maximal power of 
work-extraction [5, 6]. Naturally, t^ca is smaller than 
the Carnot upper bound 1 — Tc/Th for heat-engines. 

Heat engines have recently been studied within micro- 
scopic theories, where one is easily able to go beyond the 
quasi-equilibrium regime [7-14]. For certain classes of 
heat-engines the C A efficiency is a lower bound for the ef- 
ficiency at the maximal power of work [8-10]. This bound 
is reached at the quasi-equilibrium situation Th — > Tc in 
agreement with the finding of FTT. The result is consis- 
tent with other studies [11, 12]. 

The situation with refrigerators at a finite power is 



less clear, though [15-18]. Here maximizing the power of 
cooling does not lead to reasonable results, since there 
is an additional complementarity (not present for heat 
engines) [16]: when maximizing the power one simulta- 
neously minimizes the efficiency to zero, and vice versa. 

We study optimal regimes of finite-power refrigeration 
via a realistic model, which can be optimized over al- 
most all of its parameters. The model is quantum, but it 
admits a classical interpretation. The interest in small- 
scale refrigerators is triggered by the importance of of 
cooling processes for functioning of small devices and for 
displaying quantum features of matter [2, 14, 19-21]. 

Consider two quantum systems H and C with Hamil- 
tonians i?H and He, respectively. Each system has n 
energy levels. Initially, H and C do not interact and are 
in equilibrium at temperatures Th = 1/ I3h > Tc ^ ^/Pc'- 
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where p and cr are the initial Gibbsian density matrices 
of H and C, respectively. We write 



p = diag[r„, ...,7-1], r„ < ... < ri, 

cr = diag[s„, si], s„ < ... < si, 

Hh = diag[e„, ...,ei = 0], e„ > ... > ei. 

He = diag[/x„, ...,^1 = 0], /i„ > ... > /^i. 



(2) 
(3) 



where diag[a,..,6] is a diagonal matrix with entries 
(a,..., 6), and where without loss of generality we have 
nullified the lowest energy level of both H and C. Thus 
the overall initial density matrix is — p ® a, and the 
initial Hamiltonian Hn + He ■ 

The goal of any refrigerator is to transfer heat from the 
cooler bath to the hotter one at the expense of consuming 
work from an external source. The present refrigerator 
model functions in the following two steps; see Fig. 1. 
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FIG. 1: The refrigerator model. Two systems H and C op- 
erate between two baths at temperatures Tc < Th and are 
driven by an external potential V{t). W and Qc and Qh are, 
respectively, the work put into the overall system and the 
heats transferred from the cold bath and to the hot bath. 



1. H and C interact with each other and with external 
sources of work. The overall interaction is described via 
a time-dependent potential V{t) in the total Hamiltonian 
(?) 1 + 1 (gi He + V{t) of H -I- C. The interaction process 
is thermall isolated: V{t) is non-zero only in a short time- 
window <t < 5 and is so large there that the influence 
of all other couplings [e.g., couplings to the baths] can be 
neglected [pulsed regime] . Thus the dynamics of H -t- C 
is unitary for Q < t < 5: 



where 51; — 57(0) = p^a is the initial state defined in (1), 
Of is the final density matrix, L{ is the unitary evolution 
operator, and where T is the time-ordering operator. The 
work put into H -|- C in this process is [1] 

W = El - Ei = tr[{HH<E)l + l<E) He) (Of - O;) ], (5) 

where Ef and Ei are initial and final energies of H + C . 

2. Once the overall system H -I- C arrives at the final 
state Ofin, V{t) is switched off, and H and C (within 
some relaxation time) return back to their initial states 
(1) under influence of the hot and cold thermal baths, 
respectively. Thus the cycle is complete and can be re- 
peated again. Because the energy is conserved during 
the relaxation, the hot bath gets an amount of heat Qh, 
while the cold bath gives up the amount of heat Qc 

Qh = tr(i?H[trcOf - p]), Q, = tiiHc[a ~ trnOf]), (6) 

where trn and trc are the partial traces. Eq. (1) and the 
unitarity (4) lead to 

PhQh - PcQc = ^(Of I |0i) = tr[Of In Of - Of In 0;] , (7) 

where iS'(0f||0i) > is the relative entropy. This quan- 
tity nullifies if and only if Of = O;; otherwise it is positive. 
Eq. (7) is the Clausius inequality, with S'(0f||0i) > 
quantifying the entropy production. Eqs. (5-7) and the 
energy conservation Qh = W + Qc imply {(3c — Ph)Qc < 
f3hW, meaning that in the refrigeration regime we have 



Qc > and thus W > 0. Eq. (7) leads to the Carnot 
bound for the efficiency ^ of our refrigerator 

C = Qc/W < 0/(1 -6) = Cc, = Tc/Th < 1. (8) 

Recall that the power of refrigeration is defined as the ra- 
tio of the transferred heat to the cycle duration t, Qc/t. 
For the present model r is mainly the duration of the sec- 
ond stage, i.e., r is the relaxation time, which depends 
on the concrete physics of the system-bath coupling. For 
a weak system-bath coupling r is larger than the inter- 
nal characteristic time of H and C. In contrast, for the 
collisional system-bath interaction, t can be very short; 
see, e.g., [8] for a detailed discussion. Thus in our setup 
the cycle time r is finite, and the power of refrigeration 
Qc/t does not vanish due to a large cycle time, though 
it can vanish due to Qc — ^ 0. 

We now proceed to optimize the functioning of the 
refrigerator over the three sets of available parameters: 
the energy spacings {£k}k=2^ {Mfe}fe=2' ^^'^ ^^'^ unitary 
operators (4) [or the interaction Hamiltonian V^(i)]. 

We start by maximizing the transferred heat Qc = 
tr(77c[o' — trnOf]), which is the main characteristics of 
the refrigerator. Note that the initial energy tr[i7cO'] 
depends only on {sk}k=2- Therefore, we first choose 
{Mfe}fe=2 ^^'^ ^{t) so that the final energy tr[i7cOf] at- 
tains its minimal value equal to zero. Then we maximize 
tr[Hccr] over {efe}^^2- Note from (2, 3) 



1(g) He ^ diag[/ii , . . . , /zi, . 
Qi — pg)a — diag[siri, . . . , sir„ 



I pn : ■ ■ ■ : P"n Jj 
■ ■ 1 -^n^l 1 • ■ ■ 7 ^n^n ] 



It is clear that tr[iJcOf] = tv[HcU^liU^] goes to zero 
when, e.g., 7-2 = . . . = r„ ^ (e = £2 = . . . = £„ 00), 
while U amounts to the SWAP operation Up ® aW = 
a ® p. It is checked by a direct inspection that the maxi- 
mization of the initial energy tr[ffcf] over {efe}^^2 P'^'-'" 
duces the same structure of 71 — 1 times degenerate upper 
energy levels p, = p2 ^ ■ ■ ■ ^ P-n- Denoting 



W EE S2 = •■ = s„ = e f^^^", 



u = r2 = ..=r„ = e"'5''^(9) 



we obtain for Qc 



Qc = Tc In 



{v — u){n — 1) 



'l + {n-l)v][l + {n-l)u] 



(10) 



where according to the above discussion, Qc is maximized 
for u — > 0, and where v is to be found from maximizing 
Qcjit^o in (10) over v, i.e., v is determined via 1 -I- (n — 
l)t; + lni> = 0. Thus C can be cooled down to its ground 
state, but at a vanishing efficiency. 

For the efficiency we get for the present situation (H 
and C have n — 1 times degenerate upper levels, while U 
amounts to the SWAP operation): 



c 



^/W = 6lln[i;] (ln[u]-6lln[i>])" 



(11) 
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The maximization of Qc leads to m ^ 0, which then 
means that C in (11) goes to zero. Note that C in (11) 
reaches its maximal Carnot value 0/(1 — 9) for m = 
which nullifies the transferred heat Qc', see (10). Now 
we show that Qc tends to zero upon maximizing C, over 
all free parameters {sk}^=2^ {Mfe}fe=2 ^- Denoting 
{Nh)}^^! and {|«c)}fc=i for the eigenvectors of Hn and 
He, respectively, we note from (5, 6) that W and Qc 
feel U only via Cij\u = \{i-iijc\U\knlc)\^ ■ This matrix 
is double-stochastic [22]: Y.ijCij\ki = Y.kfii]\ki = 1- 
Conversely, for any double-stochastic matrix Cij | there 
is some unitary matrix U with matrix elements Uij^f^i, 
so that Cij\ki = |fc;P [22]. Thus, when maximizing 
various functions of W and Qc over the unitary U, we 
can directly maximize over the (n^ ~ 1)^ independent 



elements of 



n X n 



^ double stochastic matrix C. 



ij I kl ■ 



We did not find an analytic way of carrying out the 
complete maximization of C over all free parameters. 
Thus we had to rely on numerical recipes of Mathemat- 
ica 7, which for n = 1, . . . , 5 confirmed that Qc nullifies 
whenever C reaches (along any path) its maximal Carnot 
value. We believe this holds for an arbitrary n, though 
we lack any rigorous prove of this assertion. 



9uX = dyX — 0- Though we have numerically checked 
these results for n < 5 only, we again trust that they hold 
for an arbitrary n (one can, of course, always consider 
the above structure of energy spacings and U as a. useful 
ansatz). Note that u and v depend on 6* = Tc/Th- The 
efficiency C and the transferred heat Qc are given by (11) 
and (10) with u ^ u and v v; see Fig. 2. 

Since the state of H -I- C after the action of V[t) is 
(J® p, and because in the optimal regime the upper level 
for both H and C is n — 1 times degenerate, one can 
introduce non-equilibrium temperatures and T'^ for 
respectively H and C via p (x e~^'^^'^ and cr cx e^^"^'^. 
Thus, = ilni and d' = 4 In 4, where v = e-^-^ 



and 



see (9). This implies TcTh 



T'cTi As 



expected, the refrigeration condition v > u, see (10, 12), 
is equivalent to T'c < Tc < < Tl^, i.e., the cold system 
gets colder, while the hot system gets hotter. Note that 
the existence of temperatures and was not imposed, 
they emerged out of optimization. 

We eventually focus on two important limits: quasi- 
equilibrium regime 9 —^ 1, and the macroscopic regime 
Inn > 1. 

In the quasi- equilibrium regime 




FIG. 2: Solid line: efficiency of the optimized refrigerator 
versus the temperature ratio 9 = Tc/Th for n = 3; see (11). 
In the scale of this figure ({n = 2) and C,{n = 3) are almost 
indistinguishable. Dashed line: the lower bound — 1. 



Thus, neither Qc nor ^ are good target quantities for 
determining an optimal regime of refrigeration. But x = 
QcC is such a target quantity, as will be seen shortly. 
This is the most natural choice for our setup. This choice 
was also employed in [15]. Refs. [18, 19] report on other 
approaches to defining the optimal refrigeration. 

The numerical maximization of x = (Qc over {sk}^^2^ 
{pk}k=2 ^^'^ ^ been carried out for n = 1,...,5 
along the above lines. It produced the same structure: 
both H and C have n — 1 times degenerate upper levels, 
see (9), and the optimal U again corresponds to SWAP 
operation. We thus get for x = CQc [see (10, 11)] 



Xiu,v) 



Tc9{n- l){v-u) In 



2 1 



[In i-9ln i][l -I- (n - l)u][l -h (n - l)v] 



,(12) 



where u and v are found from maximizing xi^, ''^) '^ia 



X{a)\6 



Tc 



l)[l + (n-l)a]- 



In^ a, 



(13) 



maximizes for u 
dax{a)\e=i = 0: 



where a is found from 



[{n - l)a - 1] In a = 2[{n - l)a + 1] 



(14) 



We now work out the optimal u and v for 9—^1. It can 
be seen from (12) that the proper expansion parameter 
for — > 1 is a; = ^/l — 9. We represent u and v as 



fe=i 



akX 



J2^^^iak+bk-i)xK 



Substituting these expressions into duX = and dvX = 
and expanding these over x we note that and bk are 
determined by equating the 0{x'') terms: 



bo = a In — , 
a 



ai 



--mi, 

2 a' 



bi 



a 1, 
In - 24 

48 



In^ a] 



This implies for the efficiency at 9 1 {x = \/T— i 
1 _ In^a [48 + In^ a] In^ a 



c 



1 



48 



1536 



■x + Oix^). (15) 



Note that the expansion (15) does not apply for n — s- oo, 
since In a diverges in this limit; see (14). 

Eq. (15) suggests that — 1 is a lower bound for 

the efficiency at the maximal x- This is numerically 
checked to be the case for all < < 1 and all n; see 
also Fig. 2. Recalling (11) and our discussion after (12), 
we can interpret the lower bound for the efficiency as a 
lower bound on the intermediate temperature of C: 
^^^(F^ < £i < 1, i.e., cannot be too low. 
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The macroscopic regime of a n-level quantum system 
means ln(n — 1) :» 1, since for 1 weakly coupled 

particles the number of energy levels scales as . Now 
u and V in (12) arc sought via the following asymptotic 
expansions (m = n — 1) 

u = ^ pfc[mlnm]"'=, v = ^ a;fem"'=ln''TO,(16) 



where pk and Wk are found from substituting (16) into 
duX = and = and using ln(n — 1) ^ 1. In the 
first order we get pi = j^, oji 



i^l, which leads to 



c = 



i-e (1 - 6»)2 



20 ln[lnm] 
Inm 
1-6' 



ln^-__-ln 



2-6* 



+ 



Inm 



In^ m 



(17) 



Inm 



It is seen that in the macroscopic limit the efficiency con- 
verges to the Carnot value, while the transferred heat Qc 
is (in the leading order) a product of the colder temper- 
ature Tc and the "number of particles" ln(n — 1). Note 
that the obtained attainability of the Carnot bound is 
related to a finite power and a finite Qc- We see that the 
macroscopic limit does not commute with the cqiiilibrium 
limit, since the corrections in (17) diverge for ^ 1. 

Classical limit. A maximization of x = QcC can be 
carried out imposing equidistant spectra e„ = ne and 
/Lt„ = ri/U for H and C. We find that the optimal U again 
corresponds to SWAP operation. Thus, for x = x{^, v) 
we obtain 



X 



In i -61 In 4 



n{v" - m") 



{l-v){l-u) (1 - t;")(l - u") 



where v = e^""^ and il = '^'^ are found from max- 
imizing X- The efficiency ( is still given by (11). In 
the limit n 1 we get from (18): zt ^ 1, u — > 1 and 
0. Both X and C depend on one parame- 

1 



1 



{l-v"){l-u") 

ter (j) = , whose optimal value is </> = 
get in this limit: x = /, , "^/"f ^,^■> and 



c = 



Vr^. We 

1. 



Thus for a large number of equidistant energy levels 
(macro-limit) the optimal regime now implies homogene- 
ity (e 0, /i ^ 0), which is an indication of the clas- 
sical limit: under this conditional optimalization the ef- 
ficiency C is exactly equal to the [unconditional] lower 
limit CcA = - 1 



In conclusion, we have studied a model of a refriger- 
ator aiming to understand its optimal performance at a 
finite cooling power; see Fig. 1. The structure of the 
model is such that it can be optimized over almost all 
its parameters; additional constraints can and have been 
considered, though. We have confirmed an incompati- 
bility between optimizing the heat Qc transferred from 
the cold bath Tc and efficiency C: Maximizing one nul- 
lifies the other. A similar effect for a different model of 
quantum refrigerator has been reported in [20]. 



To get a balance between Qc and C wc have thus cho- 
sen to optimize their product C,Qc- This leads to a lower 
bound CcA = -^/fT§ — 1 (^ = ^) for the efficiency in 
addition to the upper Carnot bound Cc = jzrg — 1- The 
Carnot upper bound is reached (at a finite power and 
finite QcO in the macroscopic (many-level) limit of the 
model. To our knowledge such an effect has never been 
seen so far for refrigerator models. For the optimal re- 
frigerator the transferred heat Qc behaves as Qc oc Tc 
for Tc — » 0; see (10, 12, 17). This is in agreement with 
the optimal low-temperature behaviour of Qc from the 
viewpoint of the third law [21]. The lower bound CcA 
is reached in the equilibrium limit Tc ^ T^. Constrain- 
ing both systems to have homogeneous (classical) spec- 
tra, CcA is reached as an upper bound. This is just like 
within finite-time thermodynamics (FTT), when maxi- 
mizing the product of the cooling-power and efficiency 
[15], or the ratio of the efficiency and the cycle time [16]. 
In this sense Cca seems to be universal. It may play the 
same role as the Curzon-Ahlborn efficiency for heat en- 
gines 77cA, which, again, is an upper bound within FTT 
[5, 6], but appears as a lower bound for the c^nginc; models 
studied in [8-10]. Other opinions on the Curzon-Ahlborn 
efficiency for refrigerators are given in [17, 18]. 

This work has been supported by Volkswagenstiftung. 



[1] 
[2] 

[3] 

[4] 

[s: 
[6: 

[7 

[s: 

[9 

[lo; 
[11 

[12 

[13 
[14 
[15 
[16 
[17 

[18 

[19 

[20 

[21 
[22 



H.B. Callen, Thermodynamics (John Wiley, NY, 1985). 
J. Gemmer, M. Michel and G. Mahler, Quantum Ther- 
modynamics (Springer, 2004). 

V. Balzani, A. Credi and M. Venturi, Molecular Devices 
and Machines (Wiley- VCH, Weinheim, 2003). 
R.S. Berry et al. Thermodynamic Optimization of Finite- 
Time Processes (John Wiley & Sons, Chichester, 2000). 
F. Curzon and B. Ahlborn, Am. J. Phys. 43, 22 (1975). 
H. Leff, Am. J. Phys. 55, 602 (1987). 

C. Van den Broek, Phys. Rev. Lett. 95, 190602 (2005). 
R. Kosloff, J. Chem. Phys. 80, 1625 (1984). 

A.E. AUahverdyan et al, Phys. Rev. E 77, 041118 (2008). 

Z.C. Tu, J. Phys. A 41, 312003 (2008). 

Y. Izumida and K. Okuda, EPL 83, 60003 (2008). 

T. Schmiedl and U. Seifert, EPL 81, 20003 (2008). 

M. Esposito, et al, Phys. Rev. Lett. 102, 130602 (2009). 

A. E. AUahverdyan et al, J. Mod. Opt. 51, 2703 (2004). 
M.J. Henrich et al, Europhys. Lett., 76, 1057 (2006). 
Z. Yan and J. Chen, J. Phys. D 23, 136 (1990). 

S. Vclasco et al, Phys. Rev. Lett. 78, 3241 (1997). 

B. Jimenez de Cisneros et al, Phys. Rev. E 73, 057103 
(2006). 

A. Calvo Hernandez et al, Phys. Rev. E 63, 037102 
(2001). 

T. Fcldman and R. Kosloff, Phys. Rev. E, 61, 4774 
(2000). 

D. Segal, Phys. Rev. Lett. 101, 260601 (2008). 
Y. Rezek et al, EPL, 85, 30008 (2009). 

A.W. Marshall and 1. Olkin, Inequalities (Academic 
Press, New York, 1979). 



